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The ATLAS detector, built at one of the interaction points of the Large Hadron Collider, is operational and 
has been collecting data from cosmic rays. This paper describes the track based alignment of the ATLAS 
Inner Detector tracker which was performed using cosmic rays collected in 2008. The alignment algorithms 
are described and the performance of the alignment is demonstrated by showing the resulting hit residuals and 
comparing track parameters of upper and lower segments of tracks. The impact of the alignment on physics 
measurements is discussed. 



1. Introduction 

The ATLAS detector [J is a general purpose detec- 
tor built at one of the interaction points of the Large 
Hadron Collider (LHC) where proton on proton colli- 
sions with a center of mass energy of 14 TeV are ex- 
pected. The inner tracking system of ATLAS is made 
up of silicon detectors and straw drift tubes. While 
these detectors were placed with very high precision of 
the order of 100 /im, the precision required for physics 
necessitates determining the positions of the tracking 
elements to a few microns. This is only achievable by 
doing a track based alignment. The LHC has not yet 
started proton-proton collisions, however, the ATLAS 
detector is fully operational and has been collecting 
cosmic ray data. This paper describes the alignment 
achieved using this cosmic ray data and explores the 
expected impact of the alignment on the physics per- 
formance in the early data. 



2. Overview of the ATLAS Inner Detector 

The ATLAS detector consists of several systems. 
An inner tracker (Inner Detector), electromagnetic 
and hadronic calorimeters and a muon spectrometer. 
The Inner Detector is located within a solenoidal mag- 
netic field of about 2 Tesla and is made up of three 
subsystems as shown in Figure [T] 

The innermost subsystem is the pixel detector. It 
is made up of three barrel layers and three disks in 
each endcap giving at least 3 space points per track. 
There are a total of 1744 modules (1456 in the barrel 
and 144 in each endcap). The pixel cell size is 50 /im 
x 400 /mi with a corresponding resolution of 10 /im 
x 115 /im. The more precise measurement is in the 
<f> direction (in the bending direction of the magnetic 
field) and the less precise direction measures z in the 
barrel and r in the endcap. In the local frame of the 
module the directions are referred to as local x and 
local y respectively. 

The next subsystem is the SCT (Semi-Conductor 
Tracker) which consists of silicon microstrip detectors. 
In the SCT there are four barrel layers and 9 disks in 



each endcap giving 4 space points per track. It has a 
total of 4088 modules (2112 in the barrel and 988 in 
each endcap). The strip pitch is about 80 /im giving 
a resolution of 17 /im in the <j) (local x) measurement 
direction. The modules are made of two back to back 
sides which are rotated 40 mrad with respect to each 
other to give a stereo measurement. This results in a 
space point resolution of about 580 /im in z (barrel) 
or r (endcap). 

The outer subsystem is the TRT (Transition Radia- 
tion Tracker) . It is made up of straw drift tubes which 
have a diameter of 4 mm. The straws are embedded in 
a material that produces transition radiation photons 
which facilitates electron identification. On average 36 
straws are crossed per track. The barrel is segmented 
into 96 modules arranged in three rings. Each endcap 
is made up of 20 wheels, where each wheel consists of 
two four-plane structures. The resolution is 130 /im 
in the <fi measurement direction. 



3. Cosmic Ray Data Collection 

The data used to obtain the alignment presented 
here was collected from September to December 2008 
using cosmic rays. In total 2.6 million Inner Detector 
tracks were recorded with the magnetic field on and 
5 million with the magnetic field off. A number of 
these tracks, however, do not pass through all three 
subsystems. Requiring at least one hit in the SCT, 
results in 880K and 2 million tracks with field on and 
off respectively. Requiring at least one hit in the pixel 
detector, results in 190K tracks with field on and 230K 
tracks with field off. 



4. The Alignment Algorithms 

The hit residual is the distance between the track 
prediction and the measured hit. The basic approach 
to align the detector is to reduce the residuals. The 
main methods build a x 2 that is to be minimized. This 
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Figure 1: Cut-away view of the Inner Detector. 



X 2 is defined as 

X 2 = E rTV ~ lr 

tracks 

where r is a vector of the residuals and V is the cor- 
responding covariance matrix. The minimum is ob- 
tained by solving the condition: 

d X 2 /da = (2) 

where a is a vector containing the alignment con- 
stants. These are generally the 3 translations and 3 
rotations of each alignable structure (e.g. a module 
or layer). The solution to the minimization is of the 
form: 

•--(e£^|)7e£^)<*> 

\tracks / \tracks / 

= M- X b (4) 

where a full derivation can be found in [2| . The matrix 
M is a N x N matrix where N is the number of de- 
grees of freedom. If aligning all modules this will be 6 
(three translation and three rotations) times the num- 
ber of modules. In the case of the silicon detectors this 
results in a 35K x 35K matrix. The solution of this 
large system of linear equations can be obtained by 
full diagonalization (for example with LAPACK Q or 
ScaLAPACK [4]) which is computationally intensive 
or by fast solving techniques (e.g. MA27 Q) which 
can be performed on a standard workstation. The 



fast solving methods rely on the matrix being sparse 
which is generally the case. 

Solving this large matrix is referred to as the global 
X 2 approach [2, 6]. The correction of one module will 
be correlated to the movement of other modules and 
these correlations are taken into account in one go 
when solving this matrix. A few iterations are gener- 
ally required due to non linearities. 

A second approach is the local x 2 approach d, 0, § 
which ignores the correlations between modules and 
so one only needs to invert a 6 x 6 matrix for each 
module. Correlations are taken into account by iter- 
ating several times. 

The global % 2 is currently the baseline approach but 
both approaches are implemented in ATLAS and give 
consistent results. 

In addition to the x 2 minimization techniques, a ro- 
bust alignment algorithm @ has also been developed. 
This works by shifting modules according to their ob- 
served average residual offsets in an iterative fashion. 
In particular it takes advantage of the regions where 
modules overlap. 



5. The Alignment Strategy 

The alignment sequence was as follows. First the 
silicon detector (both pixel and SCT together) was 
aligned internally. A more detailed sequence of the 
silicon alignment is described below. For the pixel 
detector, survey information was available and was 
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used as a starting point for the alignment. Next the 
TRT was aligned internally. After that, the TRT was 
aligned with respect to the silicon detectors. Finally 
a "Center-of-Gravity" correction was made which ad- 
justs the overall translation and rotation of the entire 
Inner Detector such that the aligned detector has the 
same center of gravity as the nominal detector. This is 
needed as the minimization is insensitive to the overall 
translation and rotation of the Inner Detector. 

The alignment of the detector was done in a number 
of steps at different levels closely following the struc- 
tural assembly of the detector. The placement of the 
larger structures are less precisely known than that of 
the precision of the assembly of the modules within 
their substructures. The first level of alignment (re- 
ferred to as Level 1) was at the level of major subsys- 
tems which were installed as separate items. These 
are the 2 TRT endcaps and the TRT barrel, the 2 
SCT endcaps and the SCT barrel and the whole pixel 
system (however, for the silicon internal alignment the 
TRT was not included). The next level of alignment 
was at the disk and layer level. Each of the three pixel 
barrel layers were constructed from two semi-circular 
half shells. At this alignment level the 6 pixel layer 
half shells plus the 4 SCT barrel layers were aligned. 
Due to the poor illumination of the endcaps (since 
most cosmic ray muons travel predominantly in the 
vertical direction) the disks within the endcap were 
not aligned separately but rather the endcaps were 
aligned as a whole (2 SCT and 2 pixel endcaps) . 

Finally a module level alignment was done. Before 
this module level alignment was made, it was observed 
that the pixel staves have a significant bow lateral to 
the module plane which was also expected from the 
mechanical construction of the detector. To remain 
conservative not all degrees of freedom were aligned 
at the module level but rather they were restricted to 
the two degrees of freedom which were able to correct 
for the stave bow. These degrees of freedom were the 
angle about the normal of the module and the transla- 
tion in the precision measurement direction (local x) . 
At the module level only barrel modules were aligned, 
the endcaps were again kept as a whole. 



to this (local y) are shown. The figures show the dis- 
tributions before alignment where one sees wide dis- 
tributions which are not centered around zero. After 
the alignment the residual widths are reduced signif- 
icantly and well centered on zero. The widths are 
approaching that of an ideal alignment. 
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Figure 2: Unbiased residual distribution in the local x 
measurement direction of the pixel detector. 
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6. Results 
6.1. Residuals 

Since the alignment works by minimizing residuals, 
the distributions of the residuals are a key test of the 
performance of the alignment. Figures [5] - [5] show the 
residuals for the different detectors. The a quoted in 
the figures is the a of the Gaussian describing the core 
after doing a fit to a double Gaussian distribution. In 
the case of the pixel, both the high precision measure- 
ment direction (local x) and the direction orthogonal 



Figure 3: Unbiased residual distribution in the local y 
measurement direction of the pixel detector. 

If one takes the quadratic difference one estimates 
that the remaining misalignment is equivalent to a 
random displacement of the modules less than 20 /im. 
This was further tested by generating a residual mis- 
alignment set (labeled as "Day 1") with modules mis- 
placed randomly with Gaussian width of 20 /im in the 
local x and y directions which approximately repro- 
duces what is seen in data as shown in Figure [6] It is 
expected, however, that the remaining misalignment 
is not just random misplacement, but quite likely also 
includes systematic distortions. This is discussed fur- 
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Figure 4: Unbiased residual distribution in the local x 
measurement direction of the SCT detector. 
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Figure 5: Unbiased residual distribution for the TRT. 



ther in Section W?2\ 

Data from cosmic ray runs in 2009 were also re- 
cently processed using the alignments obtained with 
the 2008 data. While the residual widths were slightly 
increased the mean positions were still well centered 
on zero indicating that the detector has been quite 
stable over an extended period of time. 

6.2. Upper and Lower Track Comparison 

The reduction in width and the centering of resid- 
ual is a necessary condition to demonstrate a good 
alignment. However, it is not sufficient as a number 
of systematic distortions are insensitive or weakly con- 
strained by the minimization of the \ 2 - Cosmic ray 
tracks have the unique feature that many cross the 
upper and lower parts of the detector. One can take 
tracks that pass close to the origin (i.e. where beam 
particles collide in collision events) and split the track 
into a lower and upper segment and then refit these 
as two independent tracks. It is possible then to corn- 



Figure 6: Comparison of residual distribution in the local 
x direction for the pixel for perfect alignment, the aligned 
data and Monte-Carlo with "Day 1" residual misalignment 
of the order 20 /zm. 



pare the track parameters. This was done for tracks 
with p T > 2 GeV, \d \ < 50 mm and \z \ < 400 mm, 
where do is the transverse impact parameter measured 
with respect to the origin and zq is the z position at 
the point of closest approach to the origin. This cut 
ensures that the track has at least gone through the 
first pixel layer. 

The result of this procedure for the impact param- 
eter, where one looks at the difference of the impact 
parameter of the two tracks, is shown in Figure [7) Be- 
fore alignment there is a large shift away from zero and 
a very broad distribution. After alignment the width 
is significantly reduced. A shift of 11 /im is however 
observed, indicating there is still some improvement 
needed which is still under investigation. 
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Figure 7: Difference of transverse impact parameter of two 
tracks obtained from a cosmic ray track which is split into 
an upper and lower track and refit. 
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Since there are two tracks, an estimate of the im- 
pact parameter resolution can be obtained from the 
width of this distribution divided by \2. This re- 
sult in a resolution of 35 /zm. For comparison, the 
impact parameter resolution from collision events as 
estimated from Monte-Carlo is 20 /iin for tracks with 
p T = 5 GeV. 

Figure [5] shows the azimuthal angle, tfi, of the track 
and again good improvement is seen after alignment 
and the distribution is well centered around zero. 
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Figure 8: Difference of the phi of two tracks obtained from 
a cosmic ray track which is split into an upper and lower 
track and refit. 

Figure[5]shows the resolution oiq/p (the charge over 
momentum) as a function of the transverse momen- 
tum, px-i as obtained with this track splitting method. 
Both tracks reconstructed with the full Inner Detector 
(i.e. including the TRT) and those reconstructed with 
only the silicon detector are shown. The TRT, due to 
its large lever arm, is seen to significantly improve the 
resolution, especially at large momentum. Compar- 
ing tracks reconstructed with the full Inner Detector 
using perfectly aligned Monte-Carlo and tracks recon- 
structed in data, one can see that at low momentum 
the agreement is very good. At low momentum mul- 
tiple scattering dominates and as expected the align- 
ment is not as crucial for these tracks. However, for 
high momentum tracks there is less multiple scatter- 
ing and the impact of the remaining residual misalign- 
ment can be seen. 



7. Prospects for Collision Data Taking 
7.1. Expectations for First Collision Data 

As was mentioned above, the residual distributions 
from the cosmic ray data can be reproduced in Monte- 
Carlo by introducing random residual misalignments 
of the order of 20 fxm (see Figure ■ The random 
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Figure 9: Resolution of q/p as function of pr obtained 
from a cosmic ray track which is split into an upper and 
lower track and refit. 



Table I Size of misalignment for "Day 1" and "Day-100" 
residual misalignment sets. Modules were moved ran- 
domly with a Gaussian distribution with sigma as tab- 
ulated. 
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residual misalignment set used in simulation was re- 
ferred to as "Day 1" being an estimate of what is 
achievable on the first day of collisions given a starting 
alignment obtained solely from the cosmic ray data. 
While it is expected that we will rapidly improve the 
alignment with collision data, it is interesting to in- 
vestigate the effects that this "Day 1" alignment has 
on physics observables. Therefore some physics chan- 
nels were simulated with this misalignment set. In 
addition, a set was produced with a smaller amount 
of misalignment as an estimate of the alignment that 
might be achievable after a few months of running on 
collision data and is labeled as "Day 100" alignment. 
It should be cautioned, however, that there are many 
uncertainties to what will actually be achieved on this 
timescale. Table U shows the size of misalignment in- 
cluded in these two sets. 

The impact on the Z mass resolution from Z — ► /i/i 
events reconstructed using only the Inner Detector 
with the "Day 1" and "Day 100" alignment sets is 
shown in Figure [TOl [10fJ . As can be seen in the figure 
the impact of misalignment in the "Day 1" scenario 
shows large degradation with a contribution to the 
mass resolution (in quadrature) of 2.2 GeV. For the 
"Day 100" set there is still some degradation, however 
it is much reduced with a contribution from misalign- 
ment of about 1 GeV. 

The impact on observables of interest to B physics 
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Figure 10: Z mass resolution in Z — > fi/j, events for perfect 
alignment and for different random misalignment scenar- 
ios. 



was also investigated [To| by studying the impact 
on the mass resolution of J/tp (see Figure ITU) and 
in E>\ — > J/tpK°* events where J/ip —> fifi and 



if * — > ir^K^ . While there is some degradation in 
the "Day 1" sample, the dominance of the lower px 
tracks in these samples results in more of a contribu- 
tion from multiple scattering and less from misalign- 
ment. The effects of misalignment from the "Day 100" 
set are reduced even further to almost an insignificant 
amount. 
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Figure 11: J/tp mass resolution for perfect alignment and 
for different random alignment scenarios. 

A study of the impact of misalignment on rj-tagging 
was made with alignment sets different from those 
described above. The results of this study [ll| are 
summarized in Figure [T^l which shows the rejection of 
light jets (jets originating from gluons, up, down and 
strange quarks) for a fixed rj-tagging efficiency. The 
light jet rejection is defined as the inverse of the effi- 
ciency of a light jet being tagged as a r>jet. Different 
^-tagging algorithms were used. The IP2D algorithm 



is based on the transverse impact parameter signif- 
icance (the impact parameter divided by its error). 
The IP3D algorithm combines transverse impact pa- 
rameter significance with the longitudinal impact pa- 
rameter significance. The SV1 algorithm uses proper- 
ties of secondary vertices such as the vertex mass and 
the number of tracks used to make the vertex. The 
IP3D-SV1 combines the information from the IP3D 
and SV1 taggers. In this study four alignment sets 
were considered. For the first set, RandomlO, mod- 
ules were displaced and rotated randomly with a dis- 
tribution with width of about 10 /im. Layers and 
disks and the whole pixel layer were also displaced 
by a similar amount. Only misalignment in the pixel 
subsystem was considered. Taking into account the 
misplacements at several levels, the overall misalign- 
ment is considered to be comparable to the "Day 1" 
set above. The second set, Random5, introduces mis- 
alignment about half the size. The third set is the 
case of a perfectly aligned detector. The forth set, 
Aligned, was produced using a simulation where large 
misalignments were introduced which were typical of 
what is expected when building the detector and then 
running the actual ATLAS alignment algorithms with 
a mixture of collision-like Monte-Carlo and cosmic ray 
Monte-Carlo. For the RandomlO set the degradation 
for the IP3D-SV1 algorithm is significant (about 50%) 
but still one is able to achieve reasonable rejection 
rates. For the more realistic alignment case (Aligned) 
the degradation is only marginal with a 15% degrada- 
tion. The impact parameter based taggers appear to 
be more sensitive to misalignment than the secondary 
vertex based taggers. 
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Figure 12: Light jet rejection at a 6-tagging efficiency of 
60% for different tagging algorithms and different align- 
ment scenarios in tt events. The different fa-tagging algo- 
rithms and alignment scenarios are described in the text. 
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7.2. Systematic Distortions 



8. Conclusion 



While random misalignments cause some degrada- 
tion to the performance of physics, it is also likely 
there will be systematic distortions which could be 
difficult for the alignment to completely eliminate. 
Any deformation which will still allow helical tracks to 
be reconstructed will leave residuals unchanged. The 
minimization of the x 2 is not sensitive to such defor- 
mations. These deformations are also referred to as 
weak modes. One such deformation is a twist of the 
detector where the detector is systematically rotated 
as a function of the global z coordinate. Another such 
deformation is a curl deformation where each subse- 
quent layer is rotated progressively by larger amounts 
as a function of the radius. This particular deforma- 
tion creates a momentum bias. Tracks with one charge 
get a larger momentum and tracks with the opposite 
charge get a smaller momentum. Such a deformation 
was introduced into an alignment set and the results 
on the mass reconstruction of Z in Z — ► /i/i events was 
investigated. The size of the curl introduced was such 
that the outer most silicon layer was shifted around 
300 /im. This set is labeled as "Curl Large" . In addi- 
tion, the ATLAS alignment algorithms were run start- 
ing with a geometry with this curl misalignment. It 
was seen that the alignment procedures were able to 
remove much of the deformation indicating that the 
deformation was not a perfect weak mode. The align- 
ment set obtained after running the alignment is la- 
beled as "Curl Small" . The "Curl Large" is seen to 
give both a degradation in the resolution and intro- 
duces a bias in the mass. The "Curl Small" shows 
improvement but still some degradation with respect 
to the perfect alignment. The mass bias however is 
mostly removed. When running the alignment to pro- 
duce the 'Curl Small' set, only collision- like Monte- 
Carlo was used. The addition of events from cosmic 
rays is expected to improve the alignment further. 
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The ATLAS Inner Detector has successfully taken 
a large amount of cosmic ray data which has allowed 
a first full Inner Detector alignment to be achieved 
using the standard ATLAS alignment algorithms. 
The residuals obtained after alignment show much 
reduced widths and are well centered on zero. The 
resulting alignment is equivalent to less than a 20 /im 
residual misalignment in the barrel. The results 
from comparing track parameters from upper and 
lower segments also show good performance. The 
alignment is already at a level where it is possible 
to analyse low pr physics channels such as those of 
interest for B physics and the alignment is expected 
to rapidly improve with collision data. The tackling 
of systematic deformation is still expected to be a 
challenge although combining collision data with 
cosmic ray data is expected to help. 



References 



[i 

[2 

[3 

[4 

[5 
[6 

[9 
[10 

[11 



Figure 13: Z mass resolution for curl systematic misalign- 
ment scenarios. 



The ATLAS Collaboration, G. Aad, et al, The 

ATLAS Experiment at the CERN Large Hadron 

Collider, 2008 JINST 3 S08003. 

P. Briickman, A. Hichcur, S. Haywood, Global x 2 

approach to the Alignment of the ATLAS Silicon 

Tracking Detectors, ATL-INDET-PUB-2005-002 

(2005). 

E. Anderson, et al., LAPACK Users' Guide, 
Philadelphia PA, Society for Industrial and Ap- 
plied Mathematics, 1999. 

L.S. Blackford, et al, ScaLAPACK Users' Guide, 
Philadelphia PA, Society for Industrial and Ap- 
plied Mathematics, 1997. 

I.S. Duff and J.K. Reid, Rep. AERE R10533, 
HMSO, London, 1982. 

A. Bocci, W. Hulsbergen, TRT Alignment for 
SRI Cosmics and Beyond, ATL-INDET-PUB- 
2007-009 (2007). 

R. Hartel, diploma thesis, TU Miinchen, 2005. 
T. Gottfert, diploma thesis, Universitat 
Wiirzburg, 2006. 

F. Heinemann, Track Based Alignment of the AT- 
LAS Silicon Detectors with the Robust Alignment 
Algorithm, ATL-INDET-PUB-2007-011 (2007). 
The ATLAS Collaboration, The Impact of In- 
ner Detector Misalignments on Selected Physics, 
ATL-PHYS-PUB-2009-080 (2009). 

The ATLAS Collaboration, Expected Perfor- 
mance of the ATLAS Detector, Trigger and 
Physics, CERN-OPEN-2008-020, Geneva, 2008. 



